Astronomy & Astrophysics manuscript no. 14377 
July 28, 2010 



©ESO2010 



O 

(N 



The chemical evolution of IC10 

J. Yini'2, L. Magrini^, F. Matteucci^"*, G. A. Lanfranchi^, D. R. Gon9alves*', R. D. D. Costa'' 

Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 

Nandan Road, Shanghai, 200030, China, e-mail: jyin@shao . ac . en 

Dipartimento di Fisica, Sezione di Astronomia, Universita di Trieste, via G.B. Tiepolo 11, 34131 Trieste, Italy, e-mail: 

yin — matteucc@oats . inaf . it 

INAF - Osservatorio Astrofisico di Arcetri, Largo E. Fermi, 5, 1-50125 Firenze, Italy, e-mail: laura@arcetri . astro . it 

INAF - Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11, 3413 1 Trieste, Italy. 

NAT - Universidade Cruzeiro do Sul, R. Galvao Bueno 868, 01506-000, Sao Paulo, Brazil, e-mail: 

gustavo . lanfranchiOcruzeirodosul . edu . br 

UFRl - Observatorio do Valongo, Ladeira Pedro Antonio 43, 20080-090 Rio de laneiro, Brazil, e-mail: deniseOastro . ufr j . br 

lAG - Universidade de Sao Paulo, Rua do Matao 1226, 05508-900 Sao Paulo, Brazil, e-mail: robertoOastro . iag . usp . br 

ABSTRACT 



(N 

o 
u 

6 

u 

(N 
> 

O 
O 

in 

rn 

in 

o 
o 



X 



Context. Dwarf irregular galaxies are relatively simple unevolved objects where it is easy to test models of galactic chemical evolution. 
Aims. We attempt to determine the star formation and gas accretion history of ICIO, a local dwarf irregular for which abundance, gas, 
and mass determinations are available. 

Methods. We apply detailed chemical evolution models to predict the evolution of several chemical elements (He, O, N, S) and 
compared our predictions with the observational data. We consider additional constraints such as the present-time gas fraction, the star 
formation rate (SFR), and the total estimated mass of ICIO. We assume a dark matter halo for this galaxy and study the development of 
a galactic wind. We consider different star formation regimes: bursting and continuous. We explore different wind situations: normal 
wind, where all the gas is lost at the same rate and ii) metal-enhanced wind, where metals produced by supernovae are preferentially 
lost. We study a case without wind. We vary the star formation efficiency (SFE), the wind efficiency, and the time scale of the gas 
infall, which are the most important parameters in our models. 

Results. We find that only models with metal-enhanced galactic winds can reproduce the properties of ICIO. The star formation must 
have proceeded in bursts rather than continuously and the bursts must have been less numerous than ~10 over the whole galactic 
lifetime. Finally, ICIO must have formed by a slow process of gas accretion with a timescale of the order of 8 Gyr. 

Key words. Galaxies: abundances - Galaxies: evolution - Galaxies: irregular - Galaxies: dwarf - Galaxies: ISM - Galaxies, individual: 
ICIO 



1. Introduction 

Understanding how dwarf galaxies formed is a key aspect of 
modern observational cosmology because they are currently the 
most abundant galaxy population, and were probably even more 
abundant during the first epochs of the Universe. In addition, 
according to the hierarchical formation scenario (e.g.. White & 
Frenk [T99TT l. they are building blocks of larger structures, hence 
their study is crucial to understanding galaxy formation and evo- 
lution processes. Several questions about their formation and 
evolution still remain unsolved, among them: i) the origin of the 
differences between the irregular and spheroidal morphological 
types (e.g., Mateo 119981 ). ii) the relation between their mass and 
metallicity. Hi) the contribution of galactic and supernovae (SNe) 
winds to the evolution of the interstellar medium (ISM) metallic- 
ity, and iv) their age and star formation history (SFH) (cf., e.g., 
Gavilan et al. I2009I I. Local Universe dwarf galaxies are unique 
environments for studying these aspects. Owing to the intrin- 
sic faintness of dwarf galaxies, observational constraints to their 
evolution, such as the chemical abundances of stellar popula- 
tions of different ages, can be obtained only for nearby galaxies. 
Therefore, the large amount of information available for 
Local Group (LG) dwarf galaxies, such as star formation rate 
(SFR), H I and H2 contents, star surface density, and metallic- 
ity distribution, allows us to use them to test chemical evolu- 
tion models of different kind of galaxies (cf. Carigi et al. 120021 



Lanfranchi & Matteucci l2003l Molla & Diaz l2005T l. Elemental 
abundance ratios and their variation over cosmic time because of 
ongoing star formation (SF) are among the most important con- 
straints of chemical evolution models (e.g., Molla et al. 1996^. 
In this framework, optical spectroscopy of emission-line gas in 
galaxies is essential for deriving chemical abundances. In partic- 
ular, H II regions and planetary nebulae (PNe) provide informa- 
tion about the gas chemical composition at two different epochs 
in the galactic history: H 11 regions provide the present time 
chemical abundances, while PNe probe the distant past chem- 
istry. 

Dwarf irregular galaxies (dlrrs) are of great interest because 
they enable detailed studies of issues such as the occurrence of 
galactic winds and the chemical enrichment of the interstellar 
and intergalactic media by the direct measurement of the metal 
content of their ISM. In addition, their low level of evolution, 
as implied by their low metallicity and high gas content, makes 
these systems the most similar to primeval galaxies and, there- 
fore, the most useful for inferring the primordial galaxy condi- 
tions. It has also been proposed that dwarf irregulars could repre- 
sent the local counterparts of faint blue galaxies found in excess 
in deep galaxy counts (e.g., Ellis 1 1997b . 

In this work, we deal with the detailed chemical evolution 
of the dwarf irregular galaxy ICIO, building several chemical 
evolution models constrained by the observed properties of this 
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Table 1. Properties of ICIO 



10^ M0 (Mateo [1998]) . In TablefU the main properties of ICIO 
are summarized. 



Hubble type 


IrrlV 




1 


Distance 


740 kpc 




2 


Dynamical mass 


~ 1.7 X 10' Ms 




3 


Stellar mass 


~ (4 - 6) X 10^ 


Mo 


4,5,6 


H I mass 


~ (1.7- 2.3) X 


10** Mo 


6, 7, 8, 9, 10 


Ht mass 


1 X 10' Mo 




1 


gas fraction 


~ 0.22 - 0.36 






log(0/H) + 12 


8.2-8.3 




11, 12, 13, 14, 15 


[Fe/H]RGB 


-1.1 




16 


SFRh„ 


0.2 Mo yr-' 




7, 17 


SFRpiR 


0.05 Mo yr-i 




18, 19 


SFR1.4 GHz 


0.07 Mo yr-' 




19,20 


SFRpNc 


0.02 - 0.04 Mo 


yr-' 


14 



Refere nces. (1) van den Bergh ( [2000t ; (2) Demers et al. ([2004b : (3) 
Mateo | [T998] >: ( 4) Jarr ett et al. ( [2003) : (5) Sakai et al. ii999[ : (6) 
Vaduvescu et al. ( [20071 : (7) Glide Paz et al. ( [2003l>; (8) Huchtmeier et 
al. ( [T979ll:( 9) Wilson et a l. ([T996] >: (10) Garnett ( [2052]!: (11) Lequeux 
et al. iWJ^ : (12) Garn ett ([T990t : (13) Richer el al. (|200B : (14) MG09; 
(15 ) Lozin skaya et al. ( [2009] >: (16) Kim et al. ( [2009t; (17) Kennicutt et 
al. ([1993); (18) Melisse & Israel ( [1994) : (19) Bell ( [2003) : (20) White 
& Becker ( [1991) . 



galaxy. Our choice is motivated by recent observations (e.g., 
Sanna et al |2009l Magrini & Goncalves |2009l hereafter MG09) 
that place strong constraints on both its SFH and its chemical 
content and enrichment. 

We describe the baryonic content, the SFR, and the chemical 
abundances of ICIO in Sect. 2. The framework of the modelling 
is given in Sect. 3, whereas the results are presented in Sect. 4. 
Finally, Sect. 5 is devoted to the conclusions. 



2. Observational constraints 

ICIO is a low-mass, metal-poor, and actively star- forming 
galaxy. It was first recognized to be an extragalactic object by 
Hubble ( 1 1 936) . who defined it as 'one of the most curious ob- 
jects in the sky' . It is a member of the LG, possibly belonging 
to the Andromeda subgroup. Its location at low Galactic latitude 
implies quite uncertain determinations of its reddening (E(B-V) 
from 0.47 to 2.0) and distance (from 0.5 to 3 Mpc, see, e.g., Kim 
et aI. I2009l for a summary of distances and reddenings). Distance 
estimates locate ICIO at around 0.7 -0.8 Mpc (e.g., Demers et al. 
[20041 Kniazev et al. [20081 Sanna et al. [20081 Kim et al. [2009] l. 

2.1. Baryonic components 

The hydrogen gas associated with ICIO extends far beyond 
its optical size. The inner galaxy is located within an ex- 
tended, complex, and counter-rotating envelope (Huchtmeier 
119791 Shostak & SkilIman |I989) revealed by radio observations. 
The H I distribution has several large holes, possibly produced by 
SNe explosions (Wilcots & Miller [T998] l. 

Its atomic hydrogen mass has been estimated to be Mm ~ 
(L7-2.3)x 10*^ Mq (cf. Hu chtmei er[T979l Huchtmeier & Ritcher 
[19881 Wilson et al. [19961 GarnettHEH Vaduvescu et al. [2007] l 
and its stellar mass to be M^tar = (4 - 6) x 10^ Mq (cf. Sakai et al. 
[19991 Jarrett et al. [20031 Vaduvescu et al. [2007b . The molecular 
gas constitutes a small fraction of the total gas, having a mass 
of Mh2 ~ 1 X 10^ Mq (van den Bergh 1200(1 . Thus, its total gas 
fraction is ~ 0.22 - 0.36. Finally, its dynamical mass is 1.7 x 



2.2. The star formation rate 

Several studies of ICIO have identified stellar populations with 
different ages (Massey & Armandroff 1995^; Sakai et al. II999I 
Borissova et al. [20001 Sanna et al. [20081 ). Deep HST colour- 
magnitude diagrams have uncovered both intermediate-age, red 
clump stars, and old red horizontal branch (RGB) stars (Sanna 
et al. 12009) . Tracers of the young stellar populations are H 11 re- 
gions (Hodge & Lee II 990) and Wolf Rayet (WR) stars (e.g., 
Crowther et al. 12003b . In particular, WR stars are present in ICIO 
with the largest number per unit luminosity in the LG (Massey 
et al. |I992| Massev & Armandroff |I995) . On the other hand, the 
old/intermediate-age stellar populations are traced by the PNe 
and carbon stars (Magrini et al. 120031 Demers et al. 12004) . 

ICIO is remarkable because of its present-time high SFR. 
The current SFR was estimated from Ha, far-infrared (FIR), 
and 1.4 GHz imaging, while the intermediate age SFR was in- 
ferred from PN counting (MG09). The Ha flux implies a SFR 
of ~ 0.2 Moyr-' after an average correction for the internal 
extinction, but it could be as high as 0.6 Moyr-' if the largest 
extinction values available in the literature are adopted (Leroy 
et al. 12006) . In addition, the number of WR stars provides an 
independent estimate of the present SFR. A comparison be- 
tween the number of WRs in ICIO and a similar galaxy such as 
SMC- ~100, among them 30 (12) confirmed spectroscopically 
by Crowther et al. (2003) (Massey et al. 12003b . indicates that the 
current SFR of ICIO is 3 - 4 tim es higher than that of SMC, i.e., 
~ 0.5 Moyr-' (Leroy et al. [2006) . On the other hand, the FIR and 
radio continuum observations infer lower estimates of the SFR 
(~ 0.05 M pyr-'; Thronson et al. [T9901 White & Becker [T9921 
BeII I2003b . probably due to the escape of UV photons or/and to 
a low dust content. 

The mass of the stellar population that produced the ob- 
served PNe was estimated to be 1.2 x 10'^ Mq by MG09 from 
extended PN counting in the whole optical disk of ICIO. From 
that, they derived a lower limit to the SFR during the period of 
time when ICIO PN progenitors were born, i.e. 7-10 Gyr ago. 
The SFR value was estimated to be ~ 0.02 - 0.04 Mq yr-'. 



2.3. Ctiemical abundances 

Measuring the elemental abundances of populations with differ- 
ent ages are necessary to constrain the chemical evolution of a 
galaxy. PNe and H 11 regions are useful tools in this sense, having 
similar spectra, and thus allowing the determination with com- 
mon sets of observations and analysis techniques of the same 
elemental abundances. 

Basically, they are representative of the ISM composition at 
two different epochs in the galaxy lifetime. The progenitor stars 
of PNe, the low- and intermediate-mass stars 1 M© < M < 8 Mq, 
do not modify the composition of O, Ne, S, and Ar in the mate- 
rial ejected during and after the asymptotic giant branch (AGB) 
phase. Thus the PN abundances of these elements are charac- 
teristic of the composition of the ISM at the epoch when the 
PN progenitor was formed. On the other hand. He, N, and C 
are modified in the ejecta, because they are processed during the 
lifetime of the progenitor stars. In the case of ICIO, the observed 
PNe element abundances are consistent with those of an old pop- 
ulation. The low N/O ratios and the undetectable He n emission 
lines rule out high mass progenitor stars. Therefore, their pro- 
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genitor stars are of mass M < 1.2 Mq, and were therefore born 
during the first half of the age of the Universe (see MG09 for 
details). At the same time, the chemical abundances of H n re- 
gions are representative of the current composition of the ISM. 
Hence the abundances of the ff-elements, such as O, Ne, S, and 
Ar in H 11 regions and PNe give us a quantitative measurement of 
the ISM chemical enrichment from the epoch of PN progenitor 
formation to the present. 

Chemical abundances of H n regions in ICIO have been ob- 
tained by several authors, starting with Lequeux et al. ( |1979t . 
then recomputed by Skillman et al. ( |1989t . by Garnett et al. 
(11990b with updated atomic data and ionization correction fac- 
tors relative to those used by Lequeux et al. ( |1979t . and by the 
spectroscopic observations of Richer et al. (120011 ). The obser- 
vations by MG09 enlarged the number of studied H ii regions 
in whose spectra the electron temperature can be measured, and 
complemented them with PN observations. Spectroscopic obser- 
vations of the ionized gas were obtained by Lozinskaya et al. 
( I2009I I, who derived chemical abundances only using bright-line 
ratios. 

The chemical abundances used to constrain our chemical 
evolution model are those of MG09. We adopted this uniform 
sample to ensure consistency between the H ii region and PN 
chemical abundance determinations. Moreover, the chemical 
abundances of H ii regions observed by MG09 in common with 
previous surveys agree with their determinations. The average 
chemical abundances of PNe and H ii regions by MG09 are re- 
ported in Table |2] In the following sections, we adopt oxygen 
and sulphur abundance in PNe as tracers of the past ISM com- 
position, while we do not consider helium and nitrogen since 
their abundance might be altered by evolutionary processes oc- 
curring in the PN progenitors. We use the chemical abundances 
of H II regions to constrain our chemical evolution models at the 
present time. 



Gijnfit), and the outflow Gi,out{t), respectively. The SFR i^{t) is 
assumed to be 



m = eG(0, 



(2) 



where e is the star formation efficiency (SFE) expressed in unit 
of Gyr" ' , and is a free parameters in our work. The accretion rate 
of an element / follows the equation 



Gunf(t) = Xu„fAe-"\ 



(3) 



where Xij„f represents the abundances of the element i in the in- 
falling gas that we assume to be primordial (without metals), A is 
the normalization constant constrained to reproduce the present 
time total mass, and t is the accretion timescale. The rate of gas 
lost via galactic winds for each element is assumed to be propor- 
tional to the amount of gas present at the time t given by 



Gi.o,„{t) = WiAG(t)Xi^ou,{t), 



(4) 



where Xi^o,„{t), the abundance of the element / in the wind, is 
assumed to be identical to the abundance X,(f) in the ISM, the 
parameter A describes the efficiency of the galactic wind and is 
expressed in Gyr"' like e, w, is the efficiency weight of each 
element, hence w,/l is the effective wind efficiency of the ele- 
ment /, and A and w, are the other two free parameters in our 
model. In this work, we have studied two kinds of winds: the 
normal wind and the metal-enhanced wind. In the case of the 
normal wind, all elements are lost in the same way, i.e., w, - 1 
for all elements; however, winds, in which metals are lost pref- 
erentially relative to H and He, are referred to as "metal- en- 
hanced" winds. The assumed initial mass function (IMF) is the 
Salpeter ( 119551 ) one, and the adopted stellar yields are those of 
Woosley & Weaver d 19951 ) for massive stars and van den Hoek & 
Groenewegen ( 119971 ) for low- and intermediate-mass stars, both 
of which are metallicity-dependent. 



3. The chemical evolution model 

The chemical evolution model adopted for ICIO is the one of Yin 
et al. (12010b which is an update of the Bradamante et al. ( 11998b 
model for dwarf irregular galaxies. In this model, we can include 
different SFHs, either a sequence of short starbursts followed by 
quiescent periods or a continuous, but very mild, star formation. 
The model takes into account the feedback from SNe and stel- 
lar winds giving rise to a galactic outflow when the gas thermal 
energy exceeds the binding energy of gas. A dark matter halo is 
also considered for which the dark matter is 10 times more mas- 
sive than the luminous one, but roughly 3 times more extended 
than the luminous matter (see Yin et al. 120 1 01 for details). 

The time evolution of the fractional mass of element / in the 
gas, G/, is described by the equations 



-mXi(t) + Riit) + Gijnfit) - Gi.„,At), 



(1) 



where G/(0 = Mg(t)Xi{t) / Miitc) is the gas mass in the form 
of an element / normalized to the total baryonic mass Ml at the 
present day fc = 13 Gyr, Mg{t) is the gas mass at time t, andX,(f) 
represents the mass fraction of element / in the gas, i.e. abun- 
dance by mass. If we use the quantity G{t) = Mg(t)IMi(tc) to 
represent the total fractional mass of gas, X,(f) can be expressed 
by Gi{t)/G(t). The four items on the right hand side of Eq. ([1]) 
represent the mass change of the element / caused by the forma- 
tion of new stars iff{t)Xi{t), the material returned by either stel- 
lar winds or SNe explosions Ri{t), the infall of primordial gas 



4. Model results 

Since we do not know a priori the SFH of ICIO, we computed 
several models by varying the prescriptions for star formation. 
In particular, we developed models with either different num- 
bers of bursts or continuous SF. We refer to the models with sev- 
eral bursts of duration longer than the interburst phase as mod- 
els with "gasping SF", and models with short bursts (durations 
shorter than 1 Gyr) as "bursting" models. The galaxy is assumed 
to have been formed by the continuous infall of primordial gas. 
In addition, we tested models with and without galactic wind. 

4.1. Models without galactic winds 

For comparison purposes we developed models without feed- 
back, thus without galactic winds. In Table[3] we show the model 
parameters in the star-bursting and gasping cases, both without 
galactic winds: the SFE e, the time t corresponding to the middle 
of each starburst, the duration d of each burst, and the timescale 
for the infall t are indicated for each model. 

Models 1 and 2, assuming that ICIO has suff'ered 3 bursts (2 
bursts during the first half of the age of the Universe, 2 and 5 
Gyr, and a third occurring roughly at the present time, namely 
at 12.5 Gyr) predict acceptable present-day abundances (the left 
column of Fig. [U, if one adopts a low SFE (e = 0.1) and a 
reasonably long burst duration (c/ = 1 Gyr). For a short infall 
timescale (t = 0.5 Gyr, model 1), the abundances of different 
elements increase during the burst and maintain that level in the 
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Table 2. Average chemical abundances of PNe and H n regions in ICIO from MG09 



O/H(xl0-'') 12+log(He/H) 12+log(0/H) 12+log(N/H) 12+log(S/H) log(N/0) log(S/0) 



He/H 



PNe" 0.107±0.032 2.07±1.30 11.01±0.14 8.22+0.36 7.23±0.30 6.29±0.22 

H 11 regions" 0.107±0.031 2.31±1.59 11.02±0.12 8.28±0.29 7.41+0.39 6.76+0.31 



-1.08+0.20 -2.02+0.40 
-0.87+0.51 -1.52+0.23 



Notes. '"' Chemical abundances computed by MG09 also including upper limit electron temperature determinations. 



Table 3. Parameters of models without wind 



Model 


6 


n 


f 


d 


T 


name 


Gyr-' 




Gyr 


Gyr 


Gyr 


1 


0.1 


3 


2/5/12.5 


1/1/1 


0.5 


2 


0.1 


3 


2/5/12.5 


1/1/1 


8 


3 


0.3 


10 


9/10.5/12.95 


0.4/0. 1/.../0.1 


8 


4 


0.5 


3 


2/7/12.5 


1/7/1 


8 


5 


0.5 


10 


2/3.5/5/6.5/7.5/8.5/ 
9.5/10.5/11.5/12.75 


0.5/0.5/.../0.5 


8 



Notes. '"' the middle time of each burst. 






. Mod. 1 
Mod. 3 



5 10 

t (Gyr) 



Fig. 1. Time evolutions predicted by the models without winds. 
The left column shows the evolution of He, O, N, S abundances, 
the magenta points show the average abundances of PNe {open 
triangles, the horizontal lines from 3 to 6 Gyr denote the age 
range where PNe could be) and H ii regions {open circles) in 
Table |2l the right column shows the evolution of the total SFR, 
total gas mass, total stellar mass, total baryonic mass, and gas 
fraction, the magenta bars at 13 Gyr represent the range of the 
observational data (Table [T]i. Black-dot lines: model 1, 3 bursts 
with short infall timescale (t - 0.5 Gyr); red-dashed lines: 
model 2, 3 bursts with long infall timescale (t = 8 Gyr); green- 
solid lines: model 3, 10 bursts with long infall timescale (t = 8 
Gyr). 
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Fig. 2. Abundance ratios predicted by the models without winds. 
The cyan squares and circles are the observed abundance of PNe 
and H ii regions (MG09), and the filled magenta points are their 
average chemical abundances (Table |2|. The lines are for the 
same models as in Fig.[T] 



interburst phase, whereas for a long infall timescale (r = 8 Gyr, 
model 2) the abundances decrease during the interburst phase 
owing to the dilution of the infalling gas, thus the present-time 
abundances are lower than those predicted by the fast accretion 
model, even if all the other parameters are the same. Since there 
is no SF in these models (models 1 and 2) between 5.5 and 12.0 
Gyr and we do not predict the existence of stars with ages be- 
tween 7.5 and 1 Gyr, at variance with observations, we compiled 
another model with 10 short bursts {d ~ 0.1 Gyr) and a long 
infall timescale (t = 8 Gyr, to avoid a too high metallicity at 
the present time). This model (model 3) also shows acceptable 
abundance evolutionary tracks. 

However, all of these models form too few stars and a lot 
of gas remains, hence the predicted gas fraction (0.8 - 0.9) at 
the present time is very high and inconsistent with the observa- 
tions (0.22 - 0.36), as shown in the right column of Fig. [T] The 
abundance ratios (e.g., N/O, S/O) predicted by these models are 
consistent with the average trend of observations, but cannot ex- 
plain the scatter in the data (see Fig.|2|i. However, the spread in 
the oxygen abundance of H ii regions seems to be real, and was 
confirmed not only by the study of MG09, but also by earlier ob- 
servations (Lequeux et al. 1979, Hodge & Lee 1990 . Richer et 
al. 1200 11 1. The scale-length of metallicity inhomogeneities is of 
0.2 - 0.3 kpc, some clumps exhibit the highest measured metal- 
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Fig. 3. Time evolutions predicted by the models 4 and 5 with- 
out winds. The left column shows the evolution of He, O, N, S 
abundances, and the right column shows the evolution of the to- 
tal SFR, total gas mass, total stellar mass, total baryonic mass, 
and gas fraction. The observational constraints are the same as 
in Fig. [T] Black-dash lines: model 4, 3 bursts with higher SEE 
(e = 0.5); and red-solid lines: model 5, 10 bursts with higher 
SEE (e = 0.5). Both models adopt a long infall timescale (t = 8 
Gyr). 



licities, and an area with a more constant metallicity, around 
12H-log(0/H) = 8.2, being located between 0.2 to 0.5 kpc from 
the center (e.g. MG09). 

The origin of the inhomogeneities depends on the efficiency 
of the mixing mechanisms in dirrs. During the quiescent phases 
after the star formation episodes, the chemical elements ejected 
during the final stellar evolution phases cool until reaching the 
temperature of the ISM, and then mix with the rest of the 
ISM. The mixing is supposed to be homogeneous and efficient 
throughout the galaxy, and produced by epicyclic or radial mix- 
ing, or superbubble expansion (Roy & Kunth 1995). However, 
the mixing mechanisms operating on large scales are less ef- 
ficient for dIrrs because of the longer timescales of the gas 
mixing, and thus inhomogeneities can survive in some cases. 
In addition to ICIO, H ii regions with different chemical com- 
positions were observed in the dirr galaxies WLM (Hodge & 
Miller 1995), NGC5253 (Kobulnicky et al. 19 97), an d IC 4662 
(Hidalgo-Gamez et al. 2001). Weshow in Sect. 14.2.21 how mod- 
els that take account of chemically enhanced winds are able to 
explain the observed metallicity scatter 

If we increase the SEE to consume more gas (e - 0.5, models 
4 and 5), the system reaches very high abundances, even if the 
long infall timescale is adopted. In Eigs.[3]and |4] we show both 
the "gasping" and "bursting" model results. Both of them can 
reproduce the observed gas mass and gas fraction, but fail to 
reproduce most of the abundances. 
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Fig. 4. Abundance ratios predicted by the models 4 and 5 with- 
out winds. The points are the observed abundance of PNe and 
H II regions (MG09), same as in Eig. |2] The lines are for the 
same models as in Eig. [3] 

4.2. Models with winds 

When stellar feedback is included, galactic winds naturally de- 
velop in small galaxies. Here we describe the results of models 
with winds. As we have previously shown, if SEE is low the 
system is too gas-rich, whereas if SEE is high the system is too 
metal-rich. Therefore, the galactic wind should play an impor- 
tant role in the evolutionary history of ICIO. 

The rate of mass loss for the element / is WjAGj, where G,(f) 
is the gas mass in the form of the element / normalized to the 
total baryonic mass at the present day, A is the wind efficiency 
and is the same for all the elements, and w, is the weight of the 
element /. In Table |4] we show the parameters adopted for the 
wind models (both normal and metal-enhanced). In particular, 
wl and w2 are the models with normal winds, whereas all mod- 
els labelled mw have metal-enhanced winds. 

4.2.1. Normal winds 

In Eigs.|5]and |6] we show the results for normal winds, i.e., when 
all the elements are lost at the same rate, as opposed to metal- 
enhanced winds, i.e. when metals are preferentially lost relative 
to H and He. In other words, in normal winds we assume that w,- 
is the same for all the elements. 

In this situation, since the wind blows away a large amount 
of gas, the total mass of infalHng gas must increase to reproduce 
the observed total baryonic mass inside the galaxy at the present 
day. 

The wind efficiency is assumed to he A - 1, and both 
"bursting" and "gasping" models are tested (see Eigs.|5]and |6]l. 
According to our models, to have a final mass of 8 x 10^ Mq we 
should accrete more than 3 X lO^M©, which is 4 times higher 
than the present day total galactic mass. This implies that in our 
models an unrealistically large amount of gas is lost. 

In these models, the gas infall rates are also very high, and 
the galaxy accumulates its mass rapidly, as we can see from the 
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Table 4. Parameters of models with wind (r = 8 Gyr) 



Model 


e 


n 


f 


d 




A 


Minf 


name 


Gyr' 




Gyr 


Gyr 




Gyr-' 


10'* Mo 


wl 


0.5 


3 


2/7/12.5 


1/7/1 




1 


32 


w2 


0.5 


10 


2/3.5/5/6.5/7.5/8.5/9.5/10.5/11.5/12.75 


0.5/0.5/. 


./0.5 


1 


36 


mwl 


0.5 


3 


2/7/12.5 


1/7/1 




1 


11 


mw2 


0.5 


3 


2/7/12.5 


1/7/1 




2 


13 


mw3 


0.2 


1 


6.6 


12.8 




1.5 


15 


mw4 


0.3 


1 


8.5 


9 




1.5 


13 


mw5 


0.4 


1 


10 


6 




1.5 


12 


mw6 


0.4 


2 


5/10 


0.6/6 




1.5 


13 


mw7 


0.4 


3 


4.5/9/12.5 


3/3/1 




1.5 


13 


mw8 


0.4 


4 


4/7/10/12.5 


2/2/2/1 




1.5 


13 


mw9 


0.4 


6 


3/4.5/6/8/10/12.5 


1/1/1/1.5/1.5/1 


1.5 


13 


mwlO 


1 


4 


3/5/7/12.7 


0.8/0.8/0.8/0.6 


1.5 


13 


mwll 


1 


8 


3/4/5/6/7/9/11/12.8 


0.4/0.4/. 


./0.4 


1.5 


13 


mwl2 


1 


11 


3/4/5/6/7/8/9/10/11/12/12.9 


0.2/0.2/0.2/0.2/0.2/0.3/0.3/0.3/0.3/0.3/0.2 


1.5 


13 



Notes. '"' the middle time of each burst. 
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Fig. 5. Time evolutions predicted by the models with normal 
winds. The left column shows the evolution of He, O, N, S 
abundances, and the right column shows the evolution of the 
total SFR, total gas mass, total stellar mass, total baryonic 
mass, and gas fraction. The observational constraints are the 
same as in Fig. [1] Black-dash lines: model wl, 3 bursts (e - 
0.5, Mi„f = 32 X 10^ Mq); and red-solid lines: model w2, 10 
bursts(e - 0.5, /W,„/ = 36 x 10^ Mq). Both models adopt long 
infall timescale (r = 8 Gyr). 



evolutionary tracks of gas mass and total baryonic mass. A large 
number of stars (about one third of the present stellar mass) have 
formed during the first star formation episode resulting in a rapid 
increase in metallicity at early galactic epochs. Although the en- 
richment is decelerated by the development of the wind, the ele- 
mental abundances are much higher than those observed in PNe. 
The present gas mass and gas fraction are lower than those sug- 
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Fig. 6. Abundance ratios predicted by the models with normal 
winds. The points are the observed abundance of PNe and H n re- 
gions (MG09), same as in Fig.|2] The lines are for the same mod- 
els as in Fig.|5] 



gested by observations because of the efficient wind. As a conse- 
quence, the star formation activity decUnes rapidly and reaches 
a very low level at the present day, in contrast to the observed 
value. 

To summarize, the development of a normal wind implies 
that the galaxy should have lost an unrealistically large part of 
its accreted mass and been enriched very rapidly at early epochs. 
Therefore, it is reasonable to assume that mainly metals have 
been lost. 



4.2.2. Metal-enhanced winds 

"Metal-enhanced" winds have been suggested by several authors 
(Mac Low & Ferrara [T999l D'Ercole & Brighenti [19991 Recchi 
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Fig. 7. Time evolutions predicted by the models with metal- 
enhanced winds. The left column shows the evolution of He, O, 
N, S abundances, and the right column shows the evolution of the 
total SFR, total gas mass, total stellar mass, total baryonic mass, 
and gas fraction. The observational constraints are the same as 
in Fig.[T] Black-dot lines: model mwl, A-0, Mi„f = 8 x IO^Mq; 
red-dash lines: model mw2, A - \,Mi„f = 11 X 10** M©; and 
green- solid lines: model mw3, A - 2,M,„/ = 13 x IO^Mq. 
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Fig. 8. Abundance ratios predicted by the models with metal- 
enhanced winds. The points are the observed abundance of PNe 
and H n regions (MG09), same as in Fig.|2] The Unes are for the 
same models as in Fig.|2] 



et al. l200n Fujita et al. l2003l Recchi et al. l2008l l. Here, we as- 
sume that WH.He = 0.1, wo,s - 0.85, wn = 0.97 (for this particu- 



lar choice, we followed Romano et al. 120061) . which means that 
H and He are lost at a very low rate. 

In Figs. |2] and [8] we show the results of different wind ef- 
ficiencies A - 0, 1,2. In Fig. |7] it is clear that the stronger the 
wind, the lower the metallicity, and the deeper the minima in the 
abundance tracks occurring in interburst time. On the other hand, 
the loops in the tracks in Fig.[8]are caused by the galactic wind 
and can explain the scatter in S/O data. However, the models 
are less successful in reproducing the scatter in N/O abundance 
ratio, especially for 12-i-log(0/H)> 8.2. 

Since the observed present-time oxygen abundance, 
12+log(0/H), is mainly between 7.9 and 8.7, implying that 
the wind efficiency cannot be too high, we adopted /i = 1.5 
in this paper. This adopted value of wind efficiency implies a 
reasonable amount of infalling gas (M„,y = 13 x 10^ Mq), and 
also provides a good fit to the observed masses, as we can see in 
the next section. 



4.3. The effects of different SF regimes in tiie 
metal-eniianced-wind case 

In this section, we describe in more detail the models with metal- 
enhanced wind that we developed. We tested three types of SFH: 
continuous SF, gasping SF, and bursting SF. 

4.3.1 . Continuous star formation 

We developed and tested some models with continuous star 
formation lasting until the present time. Different SFEs were 
adopted. We found that if the SFE is low, the star formation dura- 
tion should be long enough to form the observed number of stars 
(12.8 Gyr if e = 0.2, in Figs. |9]and [TOll. On the other hand, if 
the SFE is higher, one should assume that the SF starts later and 
lasts for a shorter time (6 Gyr if e = 0.4, in Figs.|9]and [TOl i. In 
the former case, where SF is active since the very beginning of 
galaxy formation by infall, the wind begins to remove gas at ear- 
lier times (hence a higher total infall mass is required), and the 
metallicity increases very slowly, reaching an acceptable but low 
present-time value. In the latter case, since we need to reproduce 
the present time SFR, the star formation must have initiated ~6 
Gyr ago, and be very high when it begins, since the galaxy has 
already accreted a lot of gas. However, in this case no stars form 
from 7 to 10 Gyrs ago, which is inconsistent with the observa- 
tions. In Fig. [TOl we can see that the continuous star- formation 
scenario cannot explain the scatter in the metallicity data. 

4.3.2. Gasping star formation 

When we use the observed present SFR and gas fraction to con- 
strain the SFE, we find that e ~ 0.3 - 0.4, if SFRs derived 
from FIR and 1.4 GHz measurement are used. Here we adopt 
e = 0.4,M,„/ = 13x10^ Mq. 

In Figs. [TT]and[T2] we present our results for models with 
different numbers of bursts (« = 2, 3, 4, 6, where the last burst is 
ongoing at the present day). The total duration of SF is ~7 Gyr 
in all these 4 cases. 

For model mw6 (« = 2), the first short burst happens at 5 
Gyr, then star formation halts for 2 Gyr, and restarts at f = 7 Gyr. 
The first burst cannot happen too early, otherwise the minimum 
in metallicity (mainly caused by wind) will be too deep. In this 
model, we produce very few old metal-poor stars, but a bunch 
of metal-rich and younger stars, which does not reflect the star 
formation history of ICIO. 
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Fig. 9. Time evolutions predicted by the models with continu- 
ous star formation and metal-enhanced winds. The left column 
shows the evolution of He, O, N, S abundances, and the right col- 
umn shows the evolution of the total SFR, total gas mass, total 
stellar mass, total baryonic mass, and gas fraction. The obser- 
vational constraints are the same as in Fig. [T] Black-dot lines: 
model mw3, e - 0.2, d = 12.8 Gyr; red-dash lines: model 
mw4, e = 0.3, d - 9 Gyr; and green-solid lines: model mw5, 
e = 0.4,£/ = 6Gyr 



Fig. 11. Time evolutions predicted by the models with gasp- 
ing star formation and metal-enhanced winds. The left column 
shows the evolution of He, O, N, S abundances, and the right 
column shows the evolution of the total SFR, total gas mass, 
total stellar mass, total baryonic mass, and gas fraction. The ob- 
servational constraints are the same as in Fig.[T] Black-dot lines: 
model mw6, n = 2; red-dash lines: model mw7, n = 3; green- 
solid lines: model mw8, n - A; and blue-dash-dot lines: model 
mw9, n - 6 
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Fig. 10. Abundance ratios predicted by the models with contin- Fig. 12. Abundance ratios predicted by the models with gasping 

uous star formation and metal-enhanced winds. The points are star formation and metal-enhanced winds. The points are the ob- 

the observed abundance of PNe and H ii regions (MG09), and served abundance of PNe and H n regions (MG09), and are the 

are the same as in Fig.|2l The lines are for the same models as in same as in Fig.|2] The lines are for the same models as in Fig.fTT] 
Fig. 121 

In model mw7 {n - 3), there are three bursts, one during the 
first half of the galactic age, one during the second half of the 
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Fig. 13. Time evolutions predicted by the models with bursting 
star formation and metal-enhanced winds. The left column shows 
the evolution of He, O, N, S abundances, and the right column 
shows the evolution of the total SFR, total gas mass, total stellar 
mass, total baryonic mass, and gas fraction. The observational 
constraints are the same as in Fig. [T] Black-dot lines: model 
mwlO, n - A; red-dash lines: model mwll, n = 8; and green- 
solid lines: model mwl2, n - \\. 



galactic age, and the last one which remains active at the present 
time; the interburst time is about 1.5 Gyr. The results are able to 
more closely reproduce the data. 

For models mw8 (« - 4) or mw9 (n = 6), since the wind 
and the infalling gas decrease the metallicity of ISM during the 
interburst phase, the metal-poor stars can be formed if the galaxy 
experiences more star formation bursts. In these two cases, the 
models have more metal-poor stars formed in the first half of the 
galactic age and can also fit other observational data. Therefore, 
they are our best-fit models. 

In Fig. [121 we can see that the gasping star-formation sce- 
nario can reproduce the scatter in the abundance ratios, espe- 
cially in the S/O ratio, as we mentioned in Sect. 4.2.2. 



4.3.3. Bursting star formation 

In the gasping SF scenario, the model predicts a lot of stars 
whose abundances are close to the upper limit of observed PNe. 
Here a question arises: is it possible that ICIO is formed by short 
SF bursts {d < \ Gyr)? If this were true, the SFE should be 
higher than in the gasping SF scenario to form enough stars dur- 
ing these short bursts. We assume here that the SFE e - \, and 
that the total SF duration is around 3 Gyr in this case. 

In model mwlO (« — 4), three bursts occur in the first half of 
the galactic lifetime when the progenitors of PNe were born, and 
one occurs at the present time. A duration d ~ 0.8 Gyr for each 
burst is required to form enough stars. This leads to a relatively 
high metallicity after the first 3 bursts, followed by a sharp de- 
crease in the metallicity during the long interburst period before 
the last burst. In this way, extremely metal-poor stars would have 
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Fig. 14. Abundance ratios predicted by the models models with 
bursting star formation and metal-enhanced winds. The points 
are the observed abundance of PNe and H ii regions (MG09), 
and are the same as in Fig.|2] The lines are for the same models 
as in Fig. [13] 



recently formed, at variance with observations of H n regions, as 
shown in Fig. [14] To avoid these problems, we need more and 
shorter bursts. 

In model mwll (« = 8), there are five bursts before t -1 Gyr 
and three bursts after. Since the interburst time is short between 
the first five bursts, the abundances do not decrease too much. 
Lots of stars form between 3-7 Gyr and basically within the 
scatter of the data. 

In model mwl2 (« = II), we assume shorter bursts {d ~ 
0.2 - 0.3). Stars formed between 3-6 Gyr have metallicity con- 
sistent with the observed PNe, but in this model we still have 
lots of stars formed at a late age. The results of this model are 
consistent with observations. 

Compared with the former two SF regimes, the bursting one 
reproduces more closely the scatter in the abundance ratios, es- 
pecially for the N/O ratio when 10+log(O/H)>8.2, as shown in 
Fig [14] Therefore, model mwl2 can also be considered to be 
among the best ones. 

All these models have SFRs (~ 0.15 - 0.2 Mq yr"') consis- 
tent with that derived from Ha observations, but slightly higher 
than the SFR obtained from FIR and 1.4GHz measurements 
~ 0.05 - 0.07 Mq yr '(see Table [B. If we were to accept the 
SFR value derived from Ha, then the bursting SF history might 
also be acceptable for ICIO. On the other hand, the SFR val- 
ues obtained from FIR and 1 .4 GHz measurements imply that a 
gasping SF history is more likely. 



5. Conclusions 

We have studied the chemical evolution of the dwarf irregular 
galaxy ICIO for which detailed abundance data are now avail- 
able (MG09). We have adopted a chemical evolution model for 
dwarf irregulars developed by Yin et al. ( 120 1 Oi l. This model in- 
cludes detailed metallicity-dependent stellar yields for both mas- 
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sive and low and intermediate mass stars. The model also in- 
cludes feedback from SNe and stellar winds and follows the de- 
velopment of galactic winds. We have explored three cases: /) 
no winds (no feedback); //) normal winds, namely all the gas is 
lost at the same rate; and ///) metal-enhanced winds where the 
metals are lost preferentially relative to H and He. We also ex- 
plored several regimes of star formation: a) bursting mode with 
short and long bursts (a regime we called "gasping SF"); and b) 
a continuous but low SFR. 

We computed the evolutions of the abundances of several 
elements (He, N, O, S). The parameters of our model were: 1) 
the number of bursts; 2) the duration of bursts; 3) the efficiency 
of star formation; and 4) the efficiency of the galactic wind rate, 
assumed to be proportional to the amount of gas present at the 
time of the wind, which can vary from element to element in the 
case of the metal-enhanced wind. 

The observational constraints were represented by the abun- 
dances and abundance ratios of the above-mentioned elements, 
by the present-time gas mass, the SFR, and the estimated total 
mass of ICIO. 

In spite of the large number of parameters, the number of ob- 
servational constraints reproduced by our model makes us con- 
fident in concluding the following: 

- The galaxy ICIO probably formed by means of a slow gas 
accretion process with a long infall timescale of the order of 
8Gyr. 

- We do not know a priori the SFH of ICIO but our numerical 
simulations suggest that it should have experienced between 
a minimum of 4 and a maximum of 10 bursts of star forma- 
tion. If we were to trust the present-time SFR measured from 
FIR and 1 .4 GHz, the duration of the bursts should be longer 
than the interburst period, corresponding to a gasping mode 
of star formation, and the SFE should not be higher than 0.4 
Gyr ' . On the other hand, if we were to assume the present- 
time SFR inferred from Ho-, which is higher than the other 
measurements, then bursts of higher SFE (~ 1 Gyr"') but 
longer interburst periods can also be acceptable. In any case, 
we exclude a regime of continuous SF. 

- We exclude evolution without galactic winds. In this case, 
our model can reproduce neither the present-time gas frac- 
tion nor the metallicity of ICIO at the same time without in- 
cluding winds. Therefore, galactic winds must occur. 

- A metal-enhanced wind is more likely than a normal one, 
otherwise an unreasonable amount of mass is lost by the 
galaxy, making it impossible to fit the present-time amount 
of gas. In our optimal models (e.g. models mw8, mw9, and 
mwl2), a continuous wind develops following the first burst 
that carries away mostly metals, as assumed in hydrody- 
namical simulations. Both gasping and bursting models with 
metal-enhanced winds can reproduce the macroscopic prop- 
erties of ICIO and the scatter in the abundance ratios. 

Acknowledgements. J.Y. thanks the hospitahty of the Depailment of Physics 
of the University of Trieste where this work was accomphshed. J.Y., F.M. and 
L.M. acknowledge financial support from PRIN2007 from Italian Ministry of 
Reseai'ch, Prot. no. 2007JJC53X-001. J.Y. also thanks the financial support 
from the National Science Foundation of China No. 10573028, the Key Project 
No. 10833005, the Group Innovation Project No. 10821302, and 973 program No. 
2007CB815402. 



References 

Asplund, M., Grevesse, N., & Sauval, A. J. 2005, ASPC, 336, 25 
Bell, E, F. 2003, ApJ, 586, 794 



Borissova, J., Georgiev, L., Rosado, M., Kurtev, R., Bullejos, A., & Valdez- 

Gutierrez, M. 2000, A&A, 363, 130 
Bradamante, F., Matteucci, F., & D'Ercole, A. 1998, A&A 337, 338 
Carigi L., Hernandez X., & Gilmore G. 2002, MNRAS, 334, 117 
Crowther. R A.. Drissen, L.. Abbott, J. B., Royer, R, & Smartt, S. J. 2003, A&A, 

404, 483 
Demers, S., Battinelli, P, & Letarte, B. 2004, A&A 424, 125 
D'Ercole, A., & Brighenti, F, 1999. MNRAS, 309, 941 
Elhs, R. S. 1997, ARA&A, 35, 389 

Fujita, A., Martin C. L., Mac Low, M.-M., & Abel, T. 2003, ApJ, 599, 50 
Garnett, D. R. 1990, ApJ, 363, 142 
Garnett, D. R. 2002, ApJ, 581, 1019 

Gavilan, M., Molla, M., & Diaz, A. I. 2009, arXiv0903.0932 
Gil de Paz, A., Madore. B. R, & Pevunova, O. 2003, ApJS. 147, 29 
Hidalgo-Gamez, A. M., Masegosa, J., & Olofsson, K. 2001, A&A, 369, 797 
Hodge. P, & Miller., B. W. 1995, ApJ, 451, 176 
Hodge, P, & Lee, M. G. 1990, PASP 102, 26 

Hubble, E. 1936, The Realm of the Nebulae, New Haven: Yale University Press 
Huchtmeier, W. K. 1979, A&A, 75, 170 
Huchtmeier, W. K., & Ritcher, O.-G. 1988, A&A, 203, 237 
Jarrett, T. H., Chester, T., Cutri, R., Schneider, S. E., & Huchra, J. P 2003. AJ, 

125, 525 
Kennicutt, R. C, Jr, Tamblyn, P, & Congdon, C. E. 1994, ApJ, 435, 22 
Kim, M., Kim, E., Hwang, N., Lee, M. G., Im, M., Karoji, H., Noumaru, J., & 

Tanaka, I. 2009, ApJ, 703, 816 
Kniazev, A. Y, Pustilnik, S. A., & Zucker, D. B. 2008, MNRAS, 384, 1045 
Kobulnicky, H. A., Skillman, E. D., Roy, J.-R.. Walsh, J. R., & Rosa, M. R. 1997, 

ApJ, 477, 679 
Lanfranchi G. A., & Matteucci R 2003, MNRAS, 345, 71 
Lequeux, J., Peimbert, M., Rayo, J. F, Serrano, A., & Torres-Peimbert, S. 1979, 

A&A, 80, 155 
Leroy, A., Bolatto, A., Walter, R, & Blitz, L. 2006, ApJ, 643, 825 
Lozinskaya, T. A., Egorov, O. V, Moiseev, A. V., & Bizyaev, D. V. 2009, AstL, 

35, 730 
Mac Low, M.-M., & Ferrara, A. 1999, ApJ, 513, 142 
Magrini, L., Corradi, R. L. M., Greimel, R., Leisy, P., Lennon, D. J., Mampaso, 

A., Perinotto, M., Pollacco, D. L., et al. 2003, A&A, 407, 51 
Magrini, L. & Gon9alves, D. R. 2009, MNRAS, 398, 280 (MG09) 
Massey, R, Armandroft', T. E., & Conti, R S. 1992, AJ, 103, 1159 
Massey, R, & Armandroff, T. E. 1995, AJ, 109, 2470 
Massey, R, Olsen K. A. G., & Parker, J. W. 2003, PASR 1 15, 1265 
Mateo M. L. 1998, ARA&A, 36, 435 
Melisse, J. P M., & Israel, R P 1994, A&AS, 103, 391 
Molla M., & Diaz A. I., 2005, MNRAS, 358, 521 
Molla M., Ferrini R, & Diaz A. I. 1996, ApJ, 466, 668 
Recchi, S., Matteucci, F, & D'Ercole, A. 2001, MNRAS, 322, 800 
Recchi, S., Spitoni, E., Matteucci, R, & Lanfranchi, G. A. 2008, A&A, 489, 555 
Richer, M. G., Bullejos A., Borissova J., McCall, M. L., Lee, H., Kurtev, R., 

Georgiev, L., Kingsburgh, R. L., et al. 2001, A&A, 370, 34 
Romano, D., Tosi, M., & Matteucci, F 2006, MNRAS, 365, 759 
Roy, J.-R., & Kunth, D. 1995, A&A, 294, 432 
Sakai, S., Madore, B. F, & Freedman, W. L. 1999, ApJ, 511, 671 
Salpeter, E. E. 1955, ApJ, 121, 161 
Sanna, N., Bono, G., Stetson. P. B., Monelli, M., Pietrinfemi, A., Drozdovsky, 

I., Caputo, R, Cassisi, S., et al. 2008, ApJ, 688, L69 
Sanna, N.. Bono. G., Stetson, P. B., Pietrinfemi, A., Monelli, M., Cassisi, S., 

Buonanno, R., Sabbi, E., et al. 2009, ApJ, 699, L84 
Shostak, G. S., & Skillman, E. D. 1989. A&A, 214, 33 
Skillman, E. D., Kennicutt, R. C, & Hodge, P W. 1989, ApJ, 347, 875 
Thi-onson, H. A., Hunter, D. A., Casey S., & Harper, D. A. 1990, ApJ, 355, 94 
Vaduvescu, O., McCall, M. L., & Richer, M. G. 2007, AJ, 134, 604 
van den Bergh, S. 2000, The galaxies of the Local Group (Cambridge: 

Cambridge University Press) 
van den Hoek, L. B., & Groenewegen, M. A. T 1997, A&AS, 123, 305 
White, R. L., & Becker, R. H. 1992, ApJS, 79, 331 
White, S. D. M., & Frenk, C. S. 1991, ApJ, 379, 52 
Wilcots, E. M., & Miller, B. W., 1998, AJ, 116, 2363 
Wilson, C. D., Welch, D. L., Reid, I. N., Saha, A., & Hoessel, J. 1996, AL 1 1 1, 

1106 
Woosley, S. E., & Weaver, T. A., 1995, ApJS, 101, 181 
Yin, J., Matteucci, R, & Vladilo, G. 2010, submitted to A&A 



10 



